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ABSTRACT 

If galaxy clusters contain intracluster dust, the spectra of galaxies lying behind clusters should 
show attenuation by dust absorption. Wc compare the optical (3500 - 7200 A) spectra of 60,267 
luminous, early-type galaxies selected from the Sloan Digital Sky Survey to search for the signatures 
of intracluster dust in z ~0.05 clusters. We select massive, quiescent (i.e., non-star- forming) galaxies 
using an EW(Ha) < 2 A cut and consider galaxies in three bins of velocity dispersion, ranging from 
150 to 300 km s _1 . The uniformity of early-type galaxy spectra in the optical allows us to construct 
inverse- variance- weighted composite spectra with high signal-to-noise ratio (ranging from 10 2 — 10 3 ). 
We compare the composite spectra of galaxies that lie behind and adjacent to galaxy clusters and 
find no convincing evidence of dust attenuation on scales ~0.15 — 2 Mpc; we derive a generic limit 
of E{B — V) < 3 x 10 -3 mag on scales ~ 1 — 2 Mpc at the 99% confidence level, using conservative 
jackknife error bars, corresponding to a dust mass < 10 8 M@. On scales smaller than 1 Mpc this limit 
is slightly weaker, E(B -V)<8x 10 -3 mag. 



Subject headings: dust, extinction 
statistical 



1. INTRODUCTION 



galaxies: clusters: general — intergalactic medium — methods: 



Galaxy clusters are known to contain galaxies and hot 
gas, and they may contain extragalactic but intraclus- 
ter dust, or be accreting intergalactic dust from their 
neighborhoods. Indeed, there must be some intraclus- 
ter dust created by winds from the intracluster stars, 
which make up a signifi cant fraction of the t otal stel- 
lar mass in the cluster ((Ferguson et al. 19981 ). Alter- 
natively, dust could be introduced into the intracluster 
medium (ICM) thr ough such processes as cooling-flows 
flFabian et al. 1994|). galaxy or cluster mergers and col- 
lisions (Popescu et al. 2000), supernovae-driven galactic 
winds (Okaz aki et al. 1993? ) . ram pressure stripping of 
galaxies as they trav el through the intracluster medium 
Gunn fc Gott 19721) and accretion of primordial dust 



Popescu et al. 20001 ). Many of these processes have as- 
sociated timescales of order 10 8 -10 9 yr. A crucial ques- 
tion is then whether the dust thus injected into the in- 
tracluster medium can survive thermal sputtering in the 
hot gas. Typical d ust grain sputtering tim escales are 
r sp ~ 10 6 - 10 9 yr (|Draine fc Salpeter 1979T ). similar to 
the timescales of the dust-producing processes. These 
timescales imply that only the most recently injected 
dust is still surviving at any given moment in time, from 
which we conclude that the amount of dust in the ICM 
should be small and non-uniformly distributed. 

Measurements of intracluster dust have a long and 
rich history. The presence of dust was first hypoth- 
esized to explain the discrepancy between counts of 
galaxies located beh ind and adjacent to the Coma 
cluster (|Zwickv 1957|) . A first estimate of 0.4 mag 
for the magnitude of the -B-band extinction was sug- 
gested (|Zwickv 19611 1962), although infrared emis- 
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sion o f dust in the C oma Cluster was not de- 
tected (|Dwek et al. 1990l) . Zwicky's method was im- 
proved over the years as catalogs of galaxy clus- 
ters became available and B-band extinctions of or- 
der 0.2 mag were reported based on various proce- 
dures: using essentially the same approach as Zwicky, 
but using a larger sample of 15 galaxy clusters 
(|Karachentsev fc Lipovetskii 19691 ): considering color 
residuals to arrive a t the amount of absorption w ithin the 
Local Supercluster (jde Vaucouleurs et al. 197§ 1: looking 
at angular correlations among clusters and quasars, lead- 
ing to an e xtinction of about 0.12 mag over radii of sev- 
eral Mpc (Bo gart et al. 19731 ): considering correlations 
of high-rcdshift quasars with low-redshift galaxies, which 
gave evidence for dust in clusters at redshift z ~ 0.15 at a 
characteristic linear radius of bOOh^ 1 kpc, correspondin g 
to a dust sphere of mass 10 10 M Q (jBovle et al. 1988( 1. 
However, galaxy number counts are subject to a variety 
of biases and the dearth of ga laxies behind clusters co uld 
have other causes than dust (|Nollenberg et al. 20031 ). 

Correlations of quasars with nearby clusters were re- 
conside red and a -B-band ext inction of 0.15 mag was 
found (jRomani fc Maoz 1992D . A similar result was 
found to explain an excess o f higher redshift gal axies in 
nearby small galaxy groups (jGirardi et al. 19921 ). How- 
ever, comparing the color distribution of quasars be- 
hind a cluster with those in the vicinity of the clus- 
ter limited the relative reddening o f the two samples to 
E{B -V)< 0.05 mag (Ma oz 19951 ). Similar limits were 
obtained comparing color distrib utions of galaxies behind 
and removed from APM clusters (|Nollenberg et al. 20031 . 
on 1.3 Mpc sca les) and using large, elliptical galaxies 
(jFerguson 1993f) . 

There have been some, but contradictory, reports on 
dust in the central regions of clusters. A study of IRAS 
images of 56 clusters found two clusters with far-infrared 
color excesses th at could be due to 10 9 M Q of dust 
(|Wise et al. 19931 ). An average excess reddening in 10 
cooling-flow galaxy clusters of E(B — V) ~0.19 mag was 
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reported f or lines-of-sight to the center of these clusters 
(|Hu 19921 ). However, a later report found no convincing 
evidence of submillimeter dust emission in 11 cooling- 
flow clusters and set a n upper limit of 10 8 on the 
total mass of the dust (jAnnis fe Jewitt 1993T ). 

More recently, observations of six Abell clusters found 
a rough estimate of a dust mass of 10 7 M Q in the Coma 
cluster, but no evidence of dust i n the other five ob- 
served clusters dStickel et al. 20021 ); and no significant 
amount of infrared e mission from in tr acluster dust in 
Abell 2029 was found (|Bai et al. 20071 ). IChelouche et all 
(2007) reported reddening in a 0.1 < z < 0.3 sample 
of ~ 10 4 galaxy clusters by correlating the Sloan Digital 
Sky Survey cluster and quasar catalog and by comparing 
photometric and spectroscopic properties of quasars be- 
hind the clusters to those in the field. They found mean 
E(B — V) values of a few x 10~ 3 mag for sight lines pass- 
ing ~ Mpc from the clusters' centers. However, a recent 
study found no evidence of dust in 0.2 < z < 0.5 clus- 
ters from a photometric study of color excesses in several 
bands and, assuming a Galactic extinction law, derived 
an average visual ex tinction of (Ay) = 0.004±0.010 mag 
(IMuller et al. 20081) . 

In this paper we study the dust content of galaxy 
clusters by comparing the spectra of galaxies be- 
hind clusters of galaxies with those of galaxies 
not behind clusters. More specifically, we use 
the optical (3500 - 7200 A) spectra of luminous, 
early-type galaxies because they are known to 
dominate the stellar mass density o f the Universe 
(|Fukugita et al. 19981 : iHogg et al. 20021) and show great 
regularities in their proper ties (e.g.lOke fe Sandage 1968] : 

Fabcr 197i| iVisvanathan fe Sandage 1977] 

Diorgovski fe Davis 19871: iDressler et al. 1 987; 

Kormendv fe Diorgovski 19891: iBower et al. 19921 ; 

Roberts fe Havnes 1994 iBernardi et al. 2003al 2003b. 
2003c, 2003d). Their spectra show a remarkable similar- 
ity with any variation that does exist explained by th e 
environment and luminosity (Eiscnste in et al. 2 003). 
Any dust-attenuation-likc difference that can be found 
between the composite spectra of galaxies behind galaxy 
clusters and galaxies in the field can be attributed 
reliably to interactions of the galaxies' light with the 
intracluster medium. 

The difference between our study and the other pre- 
cise studies cited above, is that (1) we have more control 
over our galaxy population, (2) we are selecting on prop- 
erties that only weakly involve color and are therefore 
less likely to be biased, and (3) we have well-calibrated 
spectrophotometry of all objects. On the other hand, 
these considerations limit the size of our sample, so what 
we gain in per-object precision we lose in number of ob- 
jects, in some sense. Despite our relatively small sample, 
we obtain among the most stringent upper limits ever. 

In what follows, AB magnitudes are used throughout, 
a cosmological world model with (f2m, ^a) = (0.3, 0.7) 
is adopted, and the Hubble cons tant is H = 70 km 
s _1 Mpc -1 (|Komatsu et al. 20081) . f or the purposes of 
calculating distances (e.g. [Hogg 1999f ). 

2. DATA 

The Sloan Digital Sky Survey (SDSS) is obtaining 
u,g,r,i and z CCD imaging of 10 4 deg 2 of the northern 
Galactic sky, and from that imaging, selecting roughly 



10 6 targets for spec troscopy, most of t hem galaxies with 
r < 17.77 mag (iGunn et al. 19981: lYork et al. 20001 ; 
iStoughton et al. 2002t lAbazaiian et al. 2(H 13. 2004, 
2005). 

All the data processing, including astrometry 
(|Pier et al. 20031 ) , source identification, deblcnd- 
ing and photometry (TLupton et al. 200ll). cali- 
bration dFukugita et al. 1996t ISmith et al. 20021 ; 
llvezic et al. 20041) . spectro scopi c target selec- 
tion (lEisenstein et al. 200lt IStrauss et al. 20021 ; 
iRichards et al. 2002D. spectroscopic fiber placement 
(|Blanton et al. 2003aD . spectral data reduction and 
analysis (Schlegel & Buries 2006, in preparation; Schlegel 
2006, in preparation) are performed with automated 
SDSS software. 

We use the spectroscopic and photometric catalog from 
the NYU Value Ad ded Galaxy Catalog (NYU-VAGC; 
iBlanton et al. 20051) compiled from the SDSS Da ta Re- 
lease Four (DR4; lAdelman-McCarthv et al. 2006D . 

As we will be averaging the spectra of galaxies com- 
ing from different spectral "plates" , we depend on 
the calibration of these fluxes. The calibration pro- 
cedure is as follows (D. J. Schlegel, in preparation; 
IStoughton et al. 2002f ): Every spectral "plate" of fiber 
positions includes several faint (15.5-18.5 mag) F8 sub- 
dwarf stars. The spectrum of each standard star is 
spectrally typed by comparing with a grid of theoreti- 
cal spectra ge nerated from Kurucz model atmospheres 
(jKurucz 19921) using the spectral synthesis code SPEC- 
TRUM (|Grav fc Corballv 1994 iGrav et al. 200ll ). The 
spectra are calibrated with these F star spectra; i.e. they 
are multiplied by the function of wavelength that makes 
the F star spectra match the F star spectrophotometry 
(after correcting for Galactic reddening). This calibra- 
tion procedure produces consistent calibration at the 5% 
level 3 . In addition to this, the SDSS does not use an at- 
mospheric refraction corrector, so the effective fiber posi- 
tion on the sky shifts slightly as a function of wavelength. 
In the presence of brightness gradients, this creates a 
fluxing error. 

Redshifts are measured on the reduced spectra by an 
automated system, which models each galaxy spectrum 
as a linear combination of stellar eigenspectra (D. J. 
Schlegel, in preparation). The central velocity dispersion 
a v is determined by fitting the detailed spectral shape as 
a velocity-smoothed sum of stellar spectra (D.J. Schlegel 
& D. P. Finkbcincr, in preparation). 

The measurements of the equivalent width (EW) 
of the Ha line is m easured exactly as described in 
iQuintero et alj (|2004 ). Briefly, a linear fit of the spectral 
section to a linear combination of the mean SDSS old 
galaxy spectrum and the mean SDSS A-star spectrum 
with the locations of possible emission lines marked out 
is performed; this best fit model is then scaled down to 
have the same flux continuum as the data in the vicin- 
ity of the Ha emission line and subtracted to leave a 
continuum-subtracted line spectrum; the Ha line flux 
is then measured in a 20 A width interval centered on 
the line and converted to a rest-frame equivalent width 
with a continuum found by taking the inverse-variance- 
weighted average of two sections of the spectrum about 

3 ~ 4.4% in q-r and ~ 2.8% in r-i for galaxies, see 

http : //www . sdss . org/dr4/products/spectra/spectrophotometry . html 
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150 A in size and on either side of the emission line. This 
method fairly accurately models the absorption trough in 
the continuum, although in detail it leaves small negative 
residuals. 

The galaxy clusters used here, are 0.015 < z < 0.067 
member clusters taken from a friends-of-friends clus- 
ter catalog constructed from the SDSS DR3 main sam- 
ple g alaxies with absolu te magnitudes M .ir < -19.9 
mag (jBerlind et al. 20061) . We first consider a A/"™j n = 
10 member minimum, corresponding to a total abso- 
lute r-band magnitude M r < —21.9; for reference, the 
Virgo cluster contains 13 gal axies brighter than this limit 
(|Trentham fc Hodgkin 2 002). However, we also vary 
the minimum number of galaxies in the cluster between 
7V™j n = 5 (Mr < -20.8 mag) and A^j n = 20 (M r < -23.1 
mag). 

3. ANALYSIS 
3.1. Sample Construction 

Our technique is analogous to the classic "foreground 
screen" test of us ing stars to me asure the Galactic extinc- 
tion curve (e.g., [Calzctti 200lL and references therein). 
Here, we compare the spectra of galaxies that lie behind 
a cluster of galaxies, i.e. galaxies whose light has had to 
traverse a galaxy cluster on its way to observational as- 
tronomers, against the spectra of galaxies that do not lie 
behind a cluster. Galaxies that lie behind a galaxy clus- 
ter will constitute the target galaxies and galaxies that 
are not behind any cluster will form the control sample. 
The optical spectra of these objects must be intrinsically 
similar so that we can ascribe any measured differences to 
dust attenuation as it passes through a cluster of galax- 
ies; massive, early-t ype galaxies are know n to have ex- 
actly this property (|Eisenstein et al. 20031 ). 

Luminous, early-type galaxies are part of the red se- 
quence of galaxies and can be identified in many dif- 
ferent ways. Here, we select galaxies from the SDSS 
sample based on two properties derived from their op- 
tical spectra: the equivalent width of the Ha emission 
line, EW(Ha), and the stellar velocity dispersion, a v . 
We select quiescent (i.e., non-star-forming) galaxies us- 
ing EW(Hq)< 2 A. In addition, we restrict the sample to 
the redshift range 0.1 < z < 0.2. Figure Q] shows plots of 
EW(Hct), a v , color, and absolute magnitude of relevant 
galaxies in the SDSS sample. 

The galaxy clusters used to define our target and con- 
trol subsam ples were taken from a friends-of-friends clus- 
ter catalog (jBerlind et al. 20061 ). We emphasize that this 
catalog is not complete and that we only use part of 
the catalog to select our target and control galaxies. In 
addition to this, our galaxy spectra come from SDSS 
DR4, which has a larger coverage than SDSS DR3, out 
of which the cluster catalog was constructed. We expect 
the number of "false negatives" , i.e. galaxies that might 
be catalogued as control galaxies that are actually be- 
hind a cluster, to be small. In addition, because of the 
large number of galaxies in the control sample, the effect 
of any contaminating or misclassified galaxies will be di- 
luted in the stacking procedure. Finally, we will see that 
the principal source of errors will be due to the target 
galaxies, for which this effect is not important. 

We can now specify when a galaxy is part of the target 
subsample of galaxies and when it belongs to the control 
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Fig. 1. — Properties of galaxies in the SDSS sample: Plots of the 
quantities relevant for the sample selection. From left to right, top 
to bottom: equivalent width of the Ha line (EW(Ha)) vs. velocity 
dispersion <r„; <r„ vs. redshift; g-r color vs. EW(Ha); r-band 
absolute magnitude vs. <r„. 



subsample. For each galaxy cluster, every galaxy at a 
redshift exceeding that of the cluster that is found to be 
within 0.5 Mpc transverse distance of it, is considered to 
be behind the cluster. We define 7?-targct= 0.5 Mpc, which 
we will vary to investigate the radial dependence of any 
measured effect. Galaxies that are more than 7?. con t r oi = 
2 Mpc removed from every galaxy cluster, are classified 
as control galaxies. Any galaxy between these bounds is 
excluded from the analysis. 

Based on these criteria a primary sample with tar- 
get and control subsamples was constructed as follows: 
galaxies with EW(Ha) < 2 A, 200 < a v < 250 km s" 1 
and 0.1 < z < 0.2 were classified as target or control 
based on the values ^target = 0.5 Mpc and 7?. C ontroi = 2 
Mpc, using groups from the cluster catalog with a min- 
imum of 10 members. Figure [2] shows plots of several of 
these properties for the primary sample as well as their 
distribution on the celestial sphere. This figure illus- 
trates that the distribution of target and control galax- 
ies is approximately uniform in these properties. Several 
secondary samples were also considered, based on vari- 
ations of some of these properties, i.e. cr v , A/"™| n and 
^■target (see Table Q] for an overview of the properties of 
the various samples used). 

Added together, the total number of galaxy spectra 
that we consider in this paper is 60,267. This number 
includes the number of galaxies in the primary and sec- 
ondary samples, broken down in Table [2] by subsample 
(note, however, that there is a significant overlap of used 
spectra between certain samples), as well as the sizes of 
the target subsamples used in the determination of the 
radial dependency of dust-attenuation (see below). 
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TABLE 1 
Properties of the various samples 



Sample" 


EW(Ha) 


(T v 




\ /"miii b 
yV gal 




^-control 




(A) 


(km s — 


: ) 




(Mpc) 


(Mpc) 


Primary 


< 2 


200 < a v 


< 250 


10 


0.50 


2 


150 < a v < 200 


< 2 


150 < o-„ 


< 200 


10 


0.50 


2 


250 < cr„ < 300 


< 2 


250 < cr v 


< 300 


10 


0.50 


2 


N™= 5 


< 2 


200 < <r„ 


< 250 


5 


0.50 


2 


A£J»= 20 


< 2 


200 < cr„ 


< 250 


20 


0.50 


2 


^target = 0.25 M P C 


< 2 


200 < <T„ 


< 250 


10 


0.25 


2 


7?-targct = 1 MpC 


< 2 


200 < <t„ 


< 250 


10 


1.00 


2 



a The redshift range for all these samples is 0.1 < z < 0.2. 
b Minimum number of members to define a galaxy cluster. 

c Galaxies within a transverse distance 7?,targct of a cluster (and are at a larger 
redshift) are considered behind that cluster and make up the target subsample. 
d Galaxies more than 7^ contro i from every cluster make up the control subsample. 
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Fig. 2. — Galaxy properties for the primary sample: Plots of the 
quantities relevant for the sample selection for the primary sample 
(see Table [T] for the definition of the primary sample). The target 
subsample is represented by diamonds o, the control subsample by 
dots. From left to right, top to bottom: velocity dispersion vs. 
EW(Ha); redshift vs. EW(Ha); velocity dispersion vs. redshift; 
angular distribution of the target and control subsamples: RA vs. 
sin(DEC). 

3.2. Stacking Procedure and Comparison 

Since we are creating composite spectra from individ- 
ual spectra that are intrinsically very similar, we will 
weight the pixels of each spectrum with their inverse vari- 
ance, as this will give us the highest signal-to-noise com- 
posite spectrum. However, the overall normalization of 
each spectrum is determined by its distance from us and 
therefore the spectra are only alike when properly nor- 
malized. This can be achieved by comparing each spec- 
trum to a common template. This comparison consists 
of a simple, one-parameter fit to determine the scale fac- 
tor between the template and the spectrum to be added, 
which is used to rescale the spectrum before adding it to 
the weighted average. 

When stacking the spectra in the control subsample, 
such a template is not available a priori. Therefore, we 



have designed a procedure to self-consistently create this 
template from the constituent spectra, using a simple, 
iterative algorithm. To start off, the spectra are aver- 
aged as given, using inverse-variance-weighting, with no 
individual rescaling applied. In the next step, this first 
average features as the template which is used to scale 
each spectrum before adding it. The resultant average 
becomes the new template in the next iteration. This 
continues until the difference between successive itera- 
tions, defined as the maximum fractional difference over 
the spectrum, is less than ~ 1 part in 10 -6 . Convergence 
is achieved after about 5 iterations and can be quite dra- 
matic (see Fig. |3l which shows this convergence for the 
primary sample and the other velocity dispersion bins). 
The upshot of this iterative procedure is that the result- 
ing average is essentially the same as the template used 
to compute it, i.e. the composite control spectrum is ob- 
tained using itself as template. Finally, the composite 
spectrum is smoothed using a Gaussian filter with a ve- 
locity width equal to the maximum velocity dispersion 
of the sample. 

We use this composite control spectrum as the tem- 
plate when obtaining the composite target spectrum, 
since the spectra in these subsamples are intrinsically 
very similar. Indeed, the composite control spectrum is 
the best template to use in the sum of the target spectra, 
as it is a very high signal-to-noise representation of each 
target spectrum, and setting the overall normalization 
of each target spectrum equal to that of the compos- 
ite control spectrum will ensure that small wavelength- 
dependent, attenuation- like differences will be the pri- 
mary difference between the composite target spectrum 
and the composite control spectrum, i.e. the overall nor- 
malization of the composite spectrum of the target spec- 
tra will be very close to that of the composite control 
spectrum. The composite target spectrum obtained in 
this way is smoothed using the same Gaussian filter as 
used on the composite control spectrum. 

Finally, we fit to 
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where /target and /control is the composite spectrum of the 
target and control subsamples, r espectively, and 



is an 



arbitrary scale factor. Following ICharlot fc F all (2000), 
we parameterize the dust attenuation law as 
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Fig. 3. — Convergence of the average of the control spectra for the different velocity dispersion bins: From left to right: 200 < cr v < 250, 
150 < a v < 200, and 250 < a v < 300. 



where Ty is the y ~band optical depth. As in 
ICharlot fc FalU (|2000D . we adopt a = 0.7, which is a rea- 
sonable approximation to the shape of the Milky Way 
optical extinction curve. Note that our parameterization 
ignores the possibility that the emission of the cluster 
itself is influencing our results; however, since we expect 
the cluster light to be similar to the light of early-type 
galaxies, this should not bias our conclusions. 

3.3. Error Estimation 

The error on the individual pixels of the average spec- 
tra follows immediately from the stacking procedure. 
The inverse variance of an inverse- variance-weighted av- 
erage is given by the sum of the individual weights. The 
composite target spectra obtained have a median signal- 
to-noise ratio (SNR) of ~ 200, while the larger control 
subsamples lead to composite spectra with a median SNR 
of ~ 2500. Therefore, it is clear that the main source of 
error is due to the composite target spectrum. 

The error o~ Tv on the value of Ty is obtained by a jack- 
knife procedure. In general, a jackknife estimate of the 
variance of a statistic is obtained by dividing the sample 
into a number of (equal size) subsamples and obtaining 
the relevant statistic for each of these subsamples. The 
estimate of the variance of the statistic is then approx- 
imately equal to the variance of the values obtained for 
the subsamples (with a proportionality constant that de- 
pends on the number of subsamples, which rapidly con- 
verges to unity as the nu mber of subsamples increases; 
lEfron fc~T ibshiran i 1993f ) . Theoretically, this estimate is 
obtained from subsamples created by leaving out one of 
the "datapoints" (in our case, a "datapoint" is the spec- 
trum of a galaxy, consisting of many individual points) ; 
however, due to computational constraints, this calcula- 
tion is not always feasible. A possible way of dealing with 
this limitation is by dividing the sample into a number 
of subsets, based on a property that is unrelated to the 
relevant statistic, and creating jackk nife subsamples a s 
unions of all but one of these subsets (jShao &; Tu 1995f ). 

Whenever a full-fledged jackknife estimate was deemed 
too computationally intensive, an equal number of quan- 
tiles in declination were chosen to subdivide the sample, 
and a minimum of 200 jackknife subsamples was used 
to calculate errors in all these cases. Since the error is 



mostly due to inaccuracies in the composite target spec- 
trum, jackknife subsamples were only created using the 
target subsample; i.e. for the purpose of the jackknife 
procedure, the composite control spectrum was supposed 
to be known exactly. 

As a check on the validity of our error estimates, we 
performed a careful examination of the error estimates 
for the primary sample. Jackknife error estimates were 
obtained using different numbers of subsamples, ranging 
from ten to the maximum, 110, the number of galaxies 
in the primary target subsample, see below. Addition- 
ally, we implemented a bootstrap procedure, which works 
much in the same way as the jackknife procedure, but 
creates bootstrap samples from randomly picking "dat- 
apoints" with replacement from the set of, in this case, 
target spectra. Values up to 300 for the number of boot- 
strap samples were used. All jackknife and bootstrap 
estimates agreed on the first two significant figures of 
the error. 

3.4. Algorithm tests 

In order to examine the consistency of the stacking 
procedure and the estimation of the error on Ty , we have 
designed two algorithm-tests, which also tell us about the 
precision with which we can perform our measurement. 
First, we tested the accuracy and precision with which 
our stacking algorithm could recover a known value of 
the reddening, Ty. We selected a random subset of 200 
objects from the control subsample of the primary sam- 
ple and reddened them with a value of Ty = 0.025 with 
Gaussian noise of standard deviation 0.027 (i. e., mim- 
icking the result for the primary sample, see below). We 
found that our stacking algorithm retrieved a value of Ty 
= 0.024, with a la-error of 0.038. Similar results were 
found for different input-reddening values, however, the 
errors computed using the jackknife procedure were con- 
sistently larger, but of the same magnitude as the varia- 
tion that went in. 

To better simulate the actual parameter estimation for 
the primary sample, a larger random sample of 8000 spec- 
tra from the control subsample was chosen to be a mock- 
control-subsample and a randomly chosen subsample of 
this of one hundred spectra was artificially reddened to 
provide the mock-target-subsample. Several orders of 
magnitude of reddening were tried and we found that for 
values of Ty of ~ 10 _1 and Gaussian noise of the same 
magnitude, our algorithm returns the exact amount of 
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Fig. 4. — Result for the primary sample: average and comparison 
for the subsamples of the primary sample: the top panel shows the 
inverse-variance-weighted average spectrum of the galaxies in the 
control subsample, multiplied by the wavelength; the middle panel 
shows the same for the galaxies in the target subsample; the bottom 
panel shows a comparison of these two composite spectra (which 
should be exactly equal to zero if there were no dust absorption 
in galaxy clusters and the averaged spectra were exact), and a 
fit to a standard dust law (see the text for more details on this 
procedure). The bottom panel shows error bars (the top curve is 
an upper bound and the bottom curve a lower bound). Errors on 
the composite spectra in both the top panel and the middle panel 
are of the order of the line thickness. 

reddening and error, which means that our algorithm 
should be able to detect reddening of this magnitude 
in the real samples. Values of Ty of order 10 -2 were 
recovered by our algorithm as well, however, the jack- 
knife estimate of the error is consistently larger than the 
variation that was put in. We were unable to recover 
values of Ty of order 10 -3 and noise of the same mag- 
nitude. The computed error in this case is still of the 
order 10 -2 , which indicates that this is a lower bound 
set by measurement errors. Any errors obtained that are 
larger than this, are likely due to an intrinsic variation in 
dust absorption among clusters, while errors of this size 
and smaller arc consistent with being true measurement 
errors. 

4. RESULTS 

Figure [4] shows the result of averaging the spectra in 
the control subsample (top panel) and of averaging the 
spectra in the target subsample (middle panel) for the 
primary sample. The lower panel shows a difference plot 
of these two quantities, with the fractional difference de- 
fined as 



Qf< 



control 



Qf, 



control 



(3) 



which for the dust attenuation in equation ([T]) equals 
e~ T ( A ) — 1. A fit to equation {U gives the value of Ty, 
which together with the sizes of the samples is given in 
Table [5] Only the calculated error on the difference in 
the lower panel is shown here. This error is substantial 
and the error in Ty is, likewise, not negligible (see Table 
[2] for the error on Ty). 

Figure [5] gives the same analysis as Figure [4] for the 



TABLE 2 

Main results for the various samples 



Sample 



Ma 



Primary 


110 


21468 


0.025 


0.027 


150 < a v < 200 


122 


26531 


0.045 


0.029 


250 <a v < 300 


32 


7239 


-0.047 


0.053 


A/"™j n = 20 


425 


17993 


0.027 


0.014 


29 


22562 


-0.012 


0.059 


^target = 0.25 MpC 


31 


21468 


0.037 


0.039 


■fctar E ct= 1 MpC 


392 


21468 


-0.001 


0.013 



Note. — Results of the comparison of the average of 
the target galaxies and the average of the control galax- 
ies for the primary sample and the secondary samples. 
A/target gives the number of galaxies in the target subsam- 
ple; A/" C ontrol gives the number of galaxies in the control 
subsample; Ty is the best fit parameter to the dust law 
(cq. [TJ; ltr-errors on Ty are given in the a TV column. 

secondary samples (see also Table ^ . Inspection of this 
figure and the results in the table shows that for smaller 
target subsamples, the errors are significantly larger. 

We computed the radial dependence of Ty for the var- 
ious velocity dispersion bins by varying the value of 
T^target for these samples, while keeping the value of 
72-controi and other parameters fixed. The result is shown 
in Figure [5] The plots shown are cumulative in the sense 
that each value of 7^-targct gives dust attenuation for dust 
within a volume of radius Target- The range between 
0.15 and 1 Mpc was examined more carefully by increas- 
ing T^target in steps of 0.05 Mpc, whereas between 1 and 
2 Mpc, steps of 0.1 Mpc were used. For smaller values of 
"^•target the target subsample consists of only a few galax- 
ies, for which the composite spectrum cannot be reliably 
obtained. 

5. DISCUSSION 

The lower panel of Figure [4] shows a value of the differ- 
ence that is nearly flat over the whole wavelength range, 
and the 1-sigma error estimate on the value of Ty in 
Table [2] confirms that a null value is within the uncer- 
tainties. The other samples confirm this result: Most of 
the samples give results that are well within la of the 
null hypothesis, with a few of the samples giving a for- 
mally negative value for Ty, i.e. a negative absorption. 
None of the values of Ty are statistically significant indi- 
cators of positive or negative absorption. Regarding the 
negative absorptions, it must be remarked that the two 
most negative values are obtained from relatively small 
target subsample sizes, about 30 spectra in the target 
subsample (they both have the largest error values as 
well), which could account for an estimate that is signif- 
icantly off. The other negative value is essentially zero 
and occurs for the 7^-targct = 1 Mpc sample, which could 
simply indicate that, generically, at this radial distance 
there is no dust in galaxy clusters (see below for a dis- 
cussion of the radial dependence). 

The dust law used in equation (p} is related to the 
extinction by 

2.5 / A \-o.7 



A(X) 



Ty X 



'(55OOA) 



mag 



(4) 



In 10 V5500A/ 

which gives values of E(B — V) for the Ty values obtained 
of the order of 10~ 3 mag. Our estimate of the typically 
error on Ty, ~ 0.03, sets an upper bound on E(B — 
V) of ~ 5x 10~ 3 mag for values of 7\Ltar g ct = 0.5 Mpc. 
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Fig. 5. — Results for the secondary samples: Each of these figures shows for the various secondary samples what was shown in Figure[4] 
for the primary sample. From left to right, top to bottom: 150 < c v < 200; 250 < a v < 300; A/"™j n = 5; A/"™j n = 20; ^target = 0.25 Mpc; 

"^target = 1 Mpc. 



Comparing these errors to what we found for our mock- 
samples in Section 13.41 we see that the obtained errors 
are consistent with being true measurement errors. The 
most positive value of Ty is found for the 150 < a v < 200 
sample, which gives a reddening of E[B — V) = (8 ± 
5) x 10~ 3 mag. The most statistically significant value 
of Ty is found for the A/"™| n = 5 sample, with a value 
of Ty that is 2a from the null result, corresponding to 
E(B -V) = (5.0 ± 2.5) x 10~ 3 mag. These are still not 
very significant, but they are remarkable in that they are 
obtained for samples that have both a large target and a 
large control subsample. 

The radial dependence plots are the most instructive 
of the resulting plots as they might reveal the radial loca- 
tion of the dust content of galaxy clusters. If there were 
dust in a certain radial distance range, we would expect 
the Ty vs. target plot to be essentially zero up to the 
dust range, after which a sharp increase would occur over 
the range in which the dust occurs, followed by a grad- 
ual decline for larger values of 7^ ta rgct as more and more 
unattenuated spectra are added to the t arget subsample. 
For ex ample, to confirm the results of [Chclouch e et ail 
(|2007h . who found evidence of dust around 1 Mpc with 
no dust at smaller radii, we would expect to see a peak 
around 1 Mpc. 

In the range < 1 Mpc we do not find a consistent result 
in the three velocity dispersion bins we considered. The 
errors in this range are large because of the small num- 
ber of galaxies in the target subsamples of these bins 
(~ 100 for TZ 

target ~ 0.5 Mpc). The value of Ty in the 
range 150 < a v < 250 is positive; however, Ty is in the 



negative range in the 250 < a v < 300 bin. The amount 
of absorption rises in the interval between 7?.targct = 0.3 
Mpc and target = 0.7 Mpc for the 150 < a v < 200 bin, 
but the opposite happens for the 200 < a v < 250 bin. In 
the bottom-right panel of Figure [6l we see that the com- 
bined result of the three bins is essentially flat within the 
error range. 

The results are more consistent in the range between 1 
and 2 Mpc. All of the bins show a small amount of dust 
absorption, but all the values contain the null value in 
their error range. The overall significance of the result is 
only slightly smaller than the significance of an individual 
point, because of strong correlations between different 
T^target values (since this is a cumulative plot, the overlap 
between the different samples is significant). Combining 
the results in the 1-2 Mpc range, we find an average 
extinction E(B — V) = 0.002 mag with a significance 
of 1.5cr. At the 99% confidence level, we conclude that 
E(B - V) < 3 x 10 -3 mag. 

A similar analysis for 7^-targct ^0.5 Mpc gives an aver- 
age extinction E(B — V) = 0.004 mag, with a significance 
of 1.2cr. Therefore, in this range we can derive a limit 
E(B — V) < 8 x 10~ 3 mag. Both of these upper bounds 
are more stringent than the ones previously found. 

We can translate an upper bound on the extinction in to 
an upper bound on the dust mass using (|Kriigel 2 003) 



M dust = 1.5 x 10 fc 



E(B - V) 
3 x 10~ 3 mag V 1 Mp 



(—)' 
VI Mpc/ 



Mr. 



(5) 



which gives an upper bound Md us t ^ 10 8 M Q for ~ Mpc- 
scales. 
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FlG. 6. — Radial dependence of the dust attenuation: radial dependence of the parameter ry in the dust law (eq. [TJ for three velocity 
dispersion bins and their combined signal: top-left: 200 < cr v < 250; top-right: 150 < cr„ < 200; bottom-left: 250 < a v < 300; bottom-right: 
average of the three velocity dispersion bins. R sets 7?.target ; A/"™j n for each of these bins is 10, other parameters as for the primary sample. 
Note that the error bars in each plot are all covariant, since each datapoint is the cumulative signal from all galaxies inside 7£target- 
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The bottom-right panel of Figure [5] summarizes our 
results: over the whole range we considered, 72-targot 
~ 0.15 — 2 Mpc, we find a signal that is essentially flat 
and consistent with zero. Our upper limit on E{B — V) 
for distances between 1 and 2 Mpc from the center of a 
cluster is a few 10 -3 mag, which is consistent with the 
amount of dust extinction o bserv ed in 0.1 < z < 0.3 
clusters in iChelouche et al.l ([2007D . Future work that 
could significantly increase the size of the target subsam- 
ple could lower the upper bound found here, or confirm 
the existence of dust on the outskirts of galaxy clusters. 
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